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Abstract. We discuss the inverse problem of determining the, possibly anisotropic, conductivity of a body Q C 1™ when 
the so-called Dirichlet-to-Neumann map is locally given on a non empty portion T of the boundary d£l. We extend results 
of uniqueness and stability at the boundary, obtained by the same authors in SIAM J. Math. Anal. 33 (2001), no. 1, 153- 
171, where the Dirichlet-to-Neumann map was given on all of dQ instead. We also obtain a pointwise stability result at the 
boundary among the class of conductivities which are continuous at some point y G T. Our arguments also apply when the 
local Neumann-to-Dirichlet map is available. 

1 Introduction. 

In absence of internal sources, the electrostatic potential u in a conducting body, described by a domain 
ft C R™ , is governed by the elliptic equation 

div(crVu) = in ft, (1.1) 

where the symmetric, positive definite, matrix a — cr(x), x G ft represents the (possibly anisotropic) electric 
conductivity. The inverse conductivity problem consists of finding a when the so called Dirichlet-to-Neumann 
(D-N) map 

A ff : u\g U £ H*(dSl) — > <tVu ■ v\ au E iT^ft) 

is given for any u € if 1 (ft) solution to (1.1). Here, v denotes the unit outer normal to 9ft. If measurements 
can be taken only on one portion T of 9ft, then the relevant map is called the local Dirichlet-to-Neumann 
map. Let T be a non-empty open portion of 9ft and let us introduce the subspace of H s (9ft) 

id(r) = {/ e i^(9ft) | supp / c r}. (1.2) 

The local Dirichlet-to-Neumann map is given, in a weak formulation, as the operator such that 

(Al u, 0) = f aVu- V0, (1.3) 

1 

for any u, <p £ i? 1 (ft), u\qq, <j)\dn € Hi (T) and u is a weak solution to (1.1). 

The problem of recovering the conductivity of a body by taking measurements of voltage and current on 
its surface has came to be known as Electrical Impedance Tomography (EIT). Different materials display 
different electrical properties, so that a map of the conductivity cr(x), x G ft (ft domain in R n ) can be used 
to investigate internal properties of ft. EIT has many important applications in fields such as geophysics, 
medicine and non-destructive testing of materials. The first mathematical formulation of the inverse con- 
ductivity problem is due to A. P. Calderon [C], where he addressed the problem of whether it is possible to 
determine the (isotropic) conductivity by the D-N map. The case when measurements can be taken all over 
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the boundary has been studied extensively in the past and fundamental papers like KV1) . |KV2j , [SUj and 
[A"] show that the isotropic case can be considered solved. On the other hand the anisotropic case is still 
open and different lines of research have been pursued. One direction has been to find the conductivity up 
to a diffeomorphism which keeps the boundary fixed (see [LUj . [S], [N] and |LaU| ). 

Another direction has been the one to assume that the anisotropic conductivity is a priori known to 
depend on a restricted number of spatially-dependent parameters (see PQ, |AGj . |GL] and [L]). 
The problem of recovering the conductivity a by the knowledge of the local Dirichlet-to-Neumann map 
has been treated more recently. Lassas and Uhlmann in [LaUj recovered a connected compact real-analytic 
Riemannian manifold (M, g) with boundary by making use of the Green's function of the Laplace-Beltrami 
operator A g . See also |LaUT] . 

The procedure of reconstructing the conductivity at the boundary by local measurements has been studied 
first by Brown [B] , where the author gives a formula for reconstructing the isotropic conductivity pointwise 
at the boundary of a Lipschitz domain Q without any a priori smoothness assumption of the conductivity. 
Nakamura and Tanuma [NaTlj give a formula for the pointwise reconstruction of a conductivity continuous 
at one point x° of the boundary from the local D-N map when the boundary is C 1 near x . Under some 
additional regularity hypothesis the authors give a reconstruction formula for the normal derivatives of a 
on dft at £° € dft up to a certain order. A direct method for reconstructing the normal derivative of 
the conductivity from the local D-N map is presented in [NaT2j . The result in |NaTl] has been improved 
by Kang and Yun |KYj to an inductive reconstruction method by using only the value of a at x . The 
authors derive here also Holder stability estimates for the inverse problem to identify Riemannian metrics 
(up to isometry) on the boundary via the local D-N map. An overview on reconstructing formulas of the 
conductivity and its normal derivative can be found in |NaT3j . 

For related results of uniqueness in the interior in the case of local boundary data, we refer to Bukhgeim 
and Uhlmann [BU , Kenig, Sjostrand and Uhlmann [KSUj and Isakov [I], and, for stability, Heck and Wang 
HW . Results of stability for cases of piecewise constant conductivities and local boundary maps have also 
been obtained by Alessandrini and Vessella [AVj and by Di Cristo [D] , 

It should also be emphasized that, mainly for the applications of medical imaging, and in particular for 
breast cancer detection by EIT, rather than the local Dirichlet-to-Neumann map, one should consider the 
so-called local Neumann-to-Dirichlet (N-D) map. That is, the map associating to specified current densities 
supported on a portion T C dQ the corresponding boundary voltages, also measured on the same portion T 
of fl. 

In the present paper we study the inverse conductivity problem by local maps, concentrating on the 
issue of determining the boundary values of the conductivity and of its derivatives. We continue the line 
of investigation pursued in [AGj . by considering anisotropic unknown conductivities having the structure 
c(x) = A(x, a(x)), where A(x, t) is a known, matrix valued, function and a(x) is an unknown scalar function. 
The precise assumptions shall be illustrated in section O We improve upon the results obtained in [AGj 
under the following aspects. 

(i) The uniqueness and stability estimates are adapted to the local D-N map (Theorems 12.21 12.41 12.51 and 
Corollary 

(ii) the stability estimate at the boundary is obtained in the wider class of conductivities which are merely 
continuous in a neighborhood of some point at the boundary (Theorem I2.3|) . 

(iii) analogous results are obtained when the local D-N map is replaced by the local N-D map (Theorem 
& 

The paper is organized as follows. The main results are contained in section [5] ( subsections 12.11 12.21 for 
the local D-N, N-D maps respectively), while section [3] is devoted to the construction of singular solutions of 
equation having the same type of singularity as those in [A] but having support compactly embedded 
on a non-empty open subset of the boundary (see Theorem 13. 5|) . Proofs of the main results are given in 
section |41 (subsections 14. ll 14.21 for the local D-N, N-D maps respectively). 
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2 Main results. 

Let f2 be a domain in R n (n > 3), with Lipschitz boundary cT2. We recall, for sake of completeness, the 
definition of Lipschitz regularity of the boundary. We stick to notation already used in [AGj . 

DEFINITION 2.1. Given positive numbers L, r, h satisfying h > Lr, we say that a bounded domain G R™ 
has Lipschitz boundary if, for every x° G dil, there exists a rigid transformation of coordinates which maps 
x° into the origin, such that, setting x = (x', x n ), x' E R™ -1 , x n G R 7 we have 



for every x', y' G R™ 1 , with \x'\, \y'\ < r. 

Let us now recall the class H of functions A(x, t) introduced in [AGj . which will be considered as 
admissible conductivities. 

DEFINITION 2.2. Given p > n , the positive constants A, £, T > 0, and denoting by Sym n the class of 
n x n real valued symmetric matrices, we say that A{-, ■) € TL if the following conditions hold 



nn {x 

= {x 



(x 1 , x n ) | \x'\ < r, \x n \ < h} = 
(x', x n ) | < r, \x n \ < h, x n > f(x') }, 



where f = f(x') is a Lipschitz function defined for \x'\ < r, which satisfies 



/(0) = 
f{x')-f{y')\ < L\x'-y' 



A e W 1 ' P (Q x [A -1 , A], Symn), 
D t A G W 1,p (Cl x [A -1 , A], Sym n ), 



(2.1) 
(2.2) 




(2.3) 



A- 1 ^! 2 <A{x,t)Z'Z<\\Z 



for almost every x G Q, 

for every t£ [A -1 , A], ^ G R n . 



(2.4) 



D t A(x,t) £ ■ £ > T^l 2 , for almost every x G Q , 

for every t G [A~\ A] , £ G R". 



(2.5) 



We observe that (|2.4p is a condition of uniform ellipticity, whereas (|2.5I) is a condition of monotonicity with 
respect to the last variable t. 
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DEFINITION 2.3. For every p, < p < r we shall denote 

T p = {xe T\ dist(>, dT) > p}, (2.6) 
U p = {xeR n \dist(x,T p )<^}, (2.7) 

u* p = u p n n. (2.8) 

Here it is understood that for the empty set 0, we have dist(x, 0) = +oo. It is evident that, T being open 
and non-empty, there exists po, < po < r such that r po is also non empty. From now on we shall only 
consider values of p below po- 

Remark 2.1. We emphasize that po > is a number which depends on the choice ofTc dfl. It should be 
evident that if we choose T narrower and narrower, then po tends to and one should expect a deterioration 
in the stability estimates. 

2.1 The Dirichlet-to-Neumann map. 

We start by rigorously defining the local D-N map. We consider a given conductivity a e L°°(fl , Sym n ) 
satisfying the ellipticity condition 

A 1£ 2 < • £ < A|£| , for almost every x 6 f2, 

for every £ e M", (2.9) 

and we fix an open, non-empty subset T of dfl. We denote by < •, • > the L 2 (9f2)-pairing between H^(dil) 
and its dual H~i(dVt). 

DEFINITION 2.4. The local Dirichlet-to-Neumann map associated to a and T is the operator 

A r a : id(r) — > (ifi(r))* (2.io) 

defined by 

< Al g, -q > = / cr(x)\/u(x) ■ V(j>(x) dx, (2.11) 
J n 

i 

for any g, r\ e Hc {T), where u £ H 1 ^) is the weak solution to 

( div(cr(x)\7u(x)) = 0, in fi, 
1 u = g, on 

and 4> G H 1 ^) is any function such that 4>\dn — V i n the trace sense. 

Note that, by (|2.11|) . it is easily verified that A^ is selfadjoint. We shall denote by || • ||* the norm on the 

i i 

Banach space of bounded linear operators between H£ (T) and (Hi (T))*. 

We can now state a first stability result for the boundary values of the conductivity, assuming that the 
unknown anisotropic conductivity has the structure a(x) = A(x, a(x)) with A £ H known and a — a(x) 
unknown scalar function. 

THEOREM 2.2. (Tjipschitz stability of boundary values,). Given p > n, let Q be a bounded domain with 
Lipschitz boundary with constants L, r, h. Let T be the subset of dQ introduced above and po = Pq(T) the 
positive number introduced in Definition 2.3. If a, b are two real-valued functions satisfying 

A -1 < a[x), b(x) < A, for every (2.12) 
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II a II iv'-pp > II b || w 1 -p(Q) < ( 2 -13) 
/or some positive constant E > awe? A eH, then for any p, < p < po 

\\A(x, a(x))-A(x, b(x)) \\ L - { r p) <C\\A r A{xa) -A r A(xb) ||, . (2.14) 

-Here C > is a constant depending on n, p, L, r, h, diam(£l), p, po, X, £, J- , but not on E. 

The next Theorem improves upon the previous one, in that the regularity assumption (12.13|) is relaxed 
to mere continuity. 

THEOREM 2.3. ^Pointwise stability at the boundary,). Given p > n, let Q, T and po be as in Theorem 
\2.2l Suppose a, b are two real valued functions satisfying (12.12p and furthermore are continuous on U l p , for 
some p, < p < po. Let A G Ti., then for any x G T p 

\A(x, a(x)) - A(x, b(x))\ < C \\ A r AM - A r A{x>b) ||„, (2.15) 

where C > is a constant which depends on n, p, L, r, h, diam(fl), po, p, X, £, T only. 

Here we state our stability results for boundary values of the derivatives of the conductivity. 

THEOREM 2.4. (^Holder stability of derivatives at the boundary,). Given p, ft, V and po as in Theorem 
2.1, let a, b satisfy (2.11), (2.12) and A G Ti.. Suppose furthermore that for some p, < p < po, some 
positive integer k and some a, < a < 1 we have 

A G C k ' a (U p x [A" 1 , A] ), (2.16) 

\\ A II C k ^(U p x[\- 1 . A]) < E k- (2-17) 
\\a-b\\ ck , a{Dp) <E k . (2.18) 

Then 

\\D k (A(x,a(x))-A(x,b(x))) || LOO(fp) 

<C\\A r A ^ a) -A r A{ ^ b) ||**» (2.19) 

where 

k 

h = IT - s r-- ( 2 - 2 °) 

Here C > is a constant which depends only on n, p, L, r, h, diam(fl), p$, p, X, E, a, k, and E^. 

Under a slightly weaker assumption, we can also obtain the following uniqueness result. 

THEOREM 2.5. ('Uniqueness at the boundary,). Let p, tt, T, po, a, b, A as in Theorem 2.3. Suppose that 
for some p, < p < po and some positive integer k we have 

a-beC k (U p ). (2.21) 

If 

A r - A r 

IV A (x, a{x)) — A (x, b(x)) ' 

then 

D j (a-b)=0 on T p , for all j<k. (2.22) 

// in addition we have 

A G C k (u p x [A" 1 , A]), (2.23) 

then 

D J [A{x , a(x)) = A(x , b(xj) ) = on T p , for all j < k. (2.24) 



G 



THE LOCAL CALDERON PROBLEM 



What follows is a well-known consequence of the previous Theorem, see |KV2) and p£] for related 
arguments. 

COROLLARY 2.6. ("Uniqueness in the interiorj. Let a, b satisfy (2.11), (2.12) with p = oo. Let A G 7i 
and in addition A G TT 1,oo (0 x [A -1 , A], Symn} . Suppose that O can be partitioned into a finite number of 
Lipschitz domains, {Aj}j = % n, such that a — b is analytic on each Aj. 

If 

ly A(x,a) — ly A(x, &)> 

then we have 

A(x, a(x)) = A(x, b(x)) in O. (2.25) 

2.2 The Neumann-to-Dirichlet map. 

Let us introduce the following function spaces 



H*(dCl) = \<t> G f 0^0 

I Jan 

H-?{d£l) = {ipe H-i(dQ)\ (ib, 1) = o} 



Observe that if we consider the (global) D-N map A CT , that is the map A£ , introduced in (|2.10j) . in the special 
case when T — 90, we have that it maps onto oH 1 (50) , and, when restricted to o-ff 5 (90), it is injective 
with bounded inverse. Then we can define the global Neumann-to-Dirichlet map as follows. 

DEFINITION 2.5. The Neumann-to-Dirichlet map associated to a, N a : a H^^(dn) — > H%(dQ) is 
given by 

^ = (H**(«i>)~ 1 - ^ 

Note that N a can also be characterized as the selfadjoint operator satisfying 

< ib, N a tb >= a(x)Vu(x) ■ Vu(x) dx, (2.27) 



for every ip £ qH ? (90), where u £ i/ 1 (0) is the weak solution to the Neumann problem 

div(crVu) = 0, in O, 

aVu ■ i/\gn = ib, on 90, (2.28) 
Ian u = °- 

We are now in position to introduce the local version of the N-D map. Let T be an open portion of 90 and 
let A = SO \ T. We denote by i/ 2 (A) the closure in H 2 (90) of the space H£ (A) previously defined in 



(|1.2|) and we introduce 

H-i(T) = [ib £ H-*(dtt)\ (ib, f) = 0, for any / £ i4(A)} , (2.29) 

that is the space of distributions xb G (90) which are supported in f and have zero average on 90. The 
local N-D map is then defined as follows. 

DEFINITION 2.6. The local Neumann-to-Dirichlet map associated to a, T is the operator : o-ff ^ (T) — 
( H-i(T))* c H?(dSl) given by 

(N* i, j) = (N a 1, j), (2.30) 

for every i,j G $H _ 3(r). 
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When the local D-N map is replaced by the above defined local N-D map, completely analogous results 
to Theorems 12 .2ll2~5l and Corollary 12 . 61 could be obtained. For the sake of simplicity we state the appropriate 
version of Theorem 12.31 only. See also Remark 14.51 for further details. In what follows, we shall denote by 

|| • ||** the norm on the Banach space of bounded linear operators between qH~^(T) and (oH~i (T)j . 

THEOREM 2.7. Given p > n, let £1, T and po be as in Theorem \2.2i Suppose a, b are two real valued 
functions satisfying (|2.12p . continuous on Up, for some p, < p < po. Let A eH, then for any x E T p 

\A(x,a(x))-A(x,b(x))\<C\\N^ {xa) ~N^ xb) ||**, (2.31) 
where C > is a constant which depends on n, p, L, r, h, diam(fl), po, p, X, £, T only. 



3 Singular solutions vanishing on dft \ T. 

This section is devoted to the construction of particular solutions of equation (|l.ip . having the same type 

of singularity of those constructed in [S] but vanishing on the portion of the boundary <9f2 \ F. We consider 
the elliptic operator 

where the coefficient matrix (o~ij(x)) is symmetric and satisfies 

A _1 | £| 2 < 0fc(a:)& £j < A| £| 2 , for every x e B R , $ £ M", (3.2) 

and also 

II o~ij || w^v(b r )< E, i,j = l,...,n, (3.3) 
here p > n and A, E are positive constants. We recall the following theorem from pQ. 

THEOREM 3.1. ('Singular Solutions,). Let L satisfy (|3.1|) - (|3.3|) . For every spherical harmonic S m of 
degree m = 0, 1, 2, there exists u £ W^f(B^ \ {0}) such that 

Lu = in B R \ { 0}, (3.4) 

and furthermore, 

( Jx\ 

u{x) = log | Jx\ So r + w(x), when n = 2 and m = 0, 

\\Jx\J 

u(x) = \Jx\ 2 - n - m S m I yj^i j + w{x), otherwise, (3.5) 



where J is the positive definite symmetric matrix such that J — U (o"jj(0)) and w satisfies 

\w(x)\ + \x\\Dw(x)\<C\x\ 2 ~ n - m+a , in B R \ { }, (3.6) 



( f \D 2 w\A " < C s - n - m+a+ p , for every s, < s < R/2, 

\Js<\x\<2s ) 



(3.7) 



Here a is any number such that < a < 1 — ^, and C is a constant depending only on a, n, p, R, A, andE. 

Proof. See [1] Theorem 1.1]. ■ 
We shall also need the following. 
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LEMMA 3.2. Let the hypotheses of Theorem \3.1\ be satisfied. For every m — 0, 1, 2, . . . there exists a 
spherical harmonic S m of degree m such that the solution u given by Theorem \3.1\ also satisfies 

\Du{x)\ > \x\ for every x, < |x| < r , (3.8) 

where ro > depends only on X, E, p, m and R. 

Proof. The proof of this lemma can be obtained along the same lines as of O Lemma 3.1] and |AG( 
Section 3]. ■ 

Let us construct now solutions u of having a singularity of the same type of the above theorem in 
a point outside Q, and satisfying 

u\on G fli(r), 

in the sense of traces. To this purpose we shall make use of an augmented domain Q p . In fact, for any p, 
< p < po, one can always construct a domain f2 p with Lipschitz constants depending only on p, r, L, h 
such that 

ncn p , d^nftpCcL (3.9) 

and 

dist(x, dQ p ) > ^, for every x 6 U p . (3.10) 

If L is an operator of type (|3.1|) on fi, satisfying (|3.2|) , (|3.3|) on il, then for any p > 0, one can always 
extend the operator L to fl p in such a way so that L still satisfies (|3.2j) . (|3 . 3[) on the enlarged domain f2 p . 
As the boundary dfl is Lipschitz the unit normal vector field to the boundary may not be defined pointwise 
so we shall introduce a unitary vector field v locally defined near d£l such that: (i) v is C°° smooth, (ii) 
v is non-tangential to <9f2 (see |AG[ Section 3], for the construction procedure of the latter). The point 
z T = x° + tv, where x° G dfl, satisfies 

Ct < d(z T , DO.) < t, (3.11) 

for any r, < r < r°. Here C and t° are positive constants depending only on L, r, h \AG\ Lemma 3.3]. 
We distinguish the cases when m = or m > 0. For the case m = we recall the following asymptotic 
estimate which only requires the Holder continuity of the coefficients. 

THEOREM 3.3. Let ft and F be as in Theorem [EM For any t, < t < r 0} set z T = x° + tD, for 
some x° 6 F p and p, < p < po. If L is the operator of p.ip . with Holder continuous coefficients matrix 
u = j=i—n> with exponent < (3 < 1, t/ie Green's function G a for the Dirichlet boundary value 

problem 

J L G a {x, z T ) = — S(x — z T ), in Q p 
\ G{j (• , z T ) = 0, on dfl p 

has the form 

2-Tl 

G a (x, z T ) = C n (&et{a{z T ))y 1,2 (a-\z T ){x-z T )-(x-z T )} 2 + R(x, z T ), (3.12) 

where C n is a suitable dimensional constant and the remainder R(x, z T ) satisfies 

\R(x, z T )\ + \x- z T \ \V x R(x, z T )\ < C\x - z T \ 2 - n+a , (3.13) 

for every x G VL p , \x — z T \ < ro, where C — C(E) is a positive constant depending on E, ro is a positive 
number which depends only on the geometry of Q, and < a < [3. 

Proof of Theorem^ We refer to Mi, Chapter 1] and [MS (1.31)-(1.33)]. ■ 
As a Corollary, we also have 
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COROLLARY 3.4. The Green's function G a introduced in Theorem \3.3\ satisfies 

II G a (; z T )\\ HHn) < Ct^/ 2 , for any < r < r°, (3.14) 

where C > is a constant which only depends on diam(fl), X, L, r, h and r°. 

Proof of Corollary \S.4\ A straightforward consequence of the pointwise upper bound (|3.11[) and of the 
Caccioppoli Inequality (see, for example, [Gil Chapter 7 ]) yields 

K 

II G a \\ H i( n ) < — || G a \\ L v(p p \B c r{zr))i ( 3 - 15 ) 

where K — K(X, L, r, h, r°) is a positive constant depending only on A, L, r, h and r° and C > in (|3.15|) 
is the constant introduced in (|3 . . ■ 

For the case m > 0, we shall need stronger regularity assumptions on the coefficients. In fact, under the 
bound (HOI) we obtain 



THEOREM 3.5. Let ft and T be as in Theorem \2.%\ For any p, < p < po, let z be an arbitrary point 
in Up. For every m — 0, 1, 2,... and for every spherical harmonic S m ^ of degree m, there exists 
u E H} JU p \ {z}) n W^ C P {Q P \ {z}) such that 

Lu = in n p \{z}, (3.16) 

u = on dflp, in the trace sense (3-17) 



and it has the form 



u(x)^\J(x~z)\ 2 - n - m S m L J J^_ Z z ) )l j +v(x), (3.18) 



where J is the positive definite symmetric matrix such that J = y {o~ij(z)) r and the remainder v satisfies 
\v{x)\ + \x-z\ \Dv(x)\ <C\x- z\ 2 - n - m+a , in B p/i (z)\{z}, (3.19) 

(f \D 2 v\A" < Cs- n - m+a+ p, for every s, < 2s < p/4. (3.20) 

V J s<\x~z\<2s J 

Here a is any number such that < a < 1 — and C is a constant depending only on a, n, p, R, A, po, pandE . 

i 

Remark 3.6. Notice that, if z £ U P \Q then u G H 1 ^) and its trace satisfies u\gn £ -ff c 2 o(r). 



Proof of Theorem 13.51 With no loss of generality we can assume z = 0. Consider a positive number R 
sufficiently large so that -Br(O), the ball with centre and radius R, is such that f2 p C B R / 2 (0) CC -Br(O). 
We consider the singular solution of Theorem 13.11 on Br(0). Let us denote this solution by u m . Let wq be 
the solution to the problem 

div(crVwo) = 0, in fl p 
wq — —u m , on dftp. 

By recalling (|3.5p we get 

sup(|u m | + |Vu m |) < C u (3.21) 

consequently 

II w ||iji ( n p )< C 2 , (3.22) 
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where C± is a positive constant which depends on po, P, n and m only and C2 > depends only on po, P, n, 
m, R, L, r and ft. If we set 

u{x) = \Jx\ 2 ~ n - m S m (jJ^J + v i x ), for an y x e ( 3 - 23 ) 

and i> = it; + iuq, where w is the reminder appearing in (|3.5[) . Then m can be written as 

u(x) = u m (x) + wo(x), for any x G D, p (3.24) 

and satisfies f|3. 16[> . (|3.1T[) . moreover, by a standard interior regularity estimate 

II wo \\w*.p{B p/i {z))< C, (3-25) 

where C > depends on po, p, m and n, R, L, r and ft. Hence, recalling the bounds ()3.6|) . (|3.7|1 we obtain 
for t> = w + w , (|3~19| . (pTSOl) . ■ 

4 Proofs of the main theorems. 
4.1 The D-N map. 

The proofs of Theorem 12.41 12.51 and Corollary 12.61 follow the same line of the corresponding results in [AG] 
by replacing the singular solutions used there by those introduced in the previous Section 3 which vanish on 
dVl \ r. For this reason, we shall give the details of the proof of Theorem 12 . 2 1 only. 

Proof of Theorem \2.2\ Let x°gT p such that (a— b)(x°) —\\a—b\\ L ^^ r j and set z T — x° + tv, with < r < 
minjro, |}. Let G ai Gb be the Green's functions of Theorem 13.31 in Q p for the operators div (A(-, a(-))V-) 
and div (A(-, 6(-))V-) respectively, that is, for instance 

div(A(x, a(x))V G a (x , z T )) — —5(x — z T ), in fl p 
G a (-, z T ) = 0, on dVL p . 

and analogously for Gb- By possibly reducing r and taking < r < min{ro, |, ^}, we have that B rg (z T ) nil 
is not empty and moreover B ro (z T ) n £1 C U p . By recalling (|2.11|) and p[l (b), p. 253] we can write 



(A(x,a)- A{x,b))VG a -VG b + / (A(x, a) - A(x, b)) VG a - VG h 

= (( A A(x,a) - A A(x,b)) G a, Gb)- (4.1) 

Here and in the sequel, it is understood G a — G a (-, z T ) and analogously for Gb- By combining (14. ip with 
(f3~T2"l) and (pHS!) . we obtain 

(2 - n) 2 



(detA(a, z T )) 1/2 (det^l(6, z T )) 1/2 

f J b 2 (A(x, a) - Ape, b))Jl (x - z T ) ■ (x - z T ) 

JB T0 (z T )nn (J2(x - z T )-(x - z T )) n/2 (j^(x - Zr) ■ (x - z T )Y 

< C(E, £){ f \x-z T \^- n ^ 

\A(x,a) - A( Xl b)\\x-z T \ 2 ~ 2n 



a\B ro ( ZT ) 

A A( X , a) " A A( X , 6) II* ll Ga H H i (r) H G bH ff i (r ,' ( 4 - 2 ) 
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with J a = \J (A(z T , a(z r ))) J& = J (A(z T , b(z T ))) * and by the Holder continuity of A(ai, a(ir)), A(a;, b(x)), 
E2) yields 



C(n, £) 



(A(a; , a(a; ))- 1 - A(a; , ^(x ))- 1 ) (a: - z T ) ■ (x - z T ) 



Bro ( Zr )nn (J2(. T - Zr ).( x - Zr ))"/ 2 (j2 (a , _ Zt) . {x _ Zr) yl- 



< C{E,£)\ \x° - z T f\x - z T \ 2 ^- n) 

'B ro (z T )na 

\ x -xY\ x -z T \ 2 ^ 

B ro {z T )nn 

\a - b\\x ~ z T \ 2 - 2n 

n\B ro (z T ) 

+ II ^A( x , a )-^A { *, b ) ll*HG Q ||^ (r) ||G b ||^ (r) , (4.3) 

where C(n, £), C{E 1 £) are positive constants depending only on n, £ and on E, £ respectively. Let us 
recall 

(A(x°, aix ))- 1 - A(x°, Hx ))- 1 ) (x - z T ) ■ (x - z T ) 
D t {x Q ,t)~ x dt ) (x - z T ) ■ (x - z T ) 

la(x°) j 

-A(x°, t^Dtix , t)A(x°, t)" 1 dt) (x- z T ) ■ (x - z T ) 

la(x°) J 

> / T~ 2 \- 2 \x - z T \ 2 dt, (4.4) 



where the cllipticity and the monotonicity assumptions (12. 4p , (|2.5p had been used to obtain the lower bound 
estimate in (|4.4[) . By recalling ()4.16j) and combining (|4.3j) with (|4.4|) . we finally obtain 



C(n, £)T- 2 \- 2 \\a-b |U«. (fp) r 2 "™ < C(E, £){ T 2 - n+f3 + r 2 - n+f3 + C x ) 

+ C 2 || A^ (X; a) - A^ (X) b) ||* r 2 "™, 

where Ci is a positive constant depending only on diam(f2), A, L, r, h and r°. Consequently 

l|a-Mlz~(f,)<Cb/(r)+C 3 UA( x , a) -^A( x , b ) Ik 

where C*2 > is a constant depending only on n, A, E, £ and J 7 , C3 > is a constant depending only on 
n, A, £, J 7 , diam(£l), L, r, h and r° and /(r) — > as r — > 0. If we let r — > we obtain ()2.14|) . ■ 

We shall need three technical lemmas before we proceed with the proof of Theorem 12.31 
Given a £ L°°(f2) satisfying the ellipticity condition (|2.9p and such that it is continuous in J7 we can extend 
a to all of R™ in such a way that the ellipticity conditions are preserved and a is uniformly continuous in U p . 
We shall continue to call a such an extended function. Let us denote by ui the modulus of continuity of a in 
Up that is 

\a(x) - a(y)\ < uj(\x - y\), for any x,yeU p , (4.5) 

where u is a non negative real-valued function on R + so that u>(t) — > as t — > + . Let (j) e , e > 0, be a usual 
family of mollifying kernels with supp0 e C B £ (0). We introduce the mollification of a as 

a £ = <p £ * a{x). 
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LEMMA 4.1. For any e < p/2 we have 

\a e (x) — a(x)\ < u)(e), for any x G U p /2, (4-6) 
where a e is the mollified function of step e > 0. 
Proof of Lemma \4.1\ We have, for every x G R™, 

a e (x)-a(x)= <j> e (x - y) (a(y) - a(x)) dy, (4.7) 

J \y — x\<e 

and when x G U p / 2 , £ < p/2, \y — x\ < e implies y G U p , hence 



\a e {x) - a{x)\ < / \<j) e (x -y)\u(e) dy 

J\y-x\<s 

= lu(e), for any x G U p / 2 - 



LEMMA 4.2. With the same assumptions as above 

V A(-,a e (-)) - A A(-,a(-)) H* 

1 

Proof of Lemma \J7^ Let <\> € 77 c 2 (r), < e < p/2 and take u, u e £ 1 (SI) solutions to the problems 



div(,4(:r, a)Vu) =0 in 1! 

m = </> on 917 

and 

div(A(a;, a £ )Vu e ) = in Q 



u e = cf> on dQ. 

respectively, then by p. lip (see jS, (b), p. 253]) we have 



(( A A(,a e (.)) " A A(-M-)))^ $ = IJ A ( X > fl£ ) " fl )) V ^ • VU 

(A(x, a e ) — A(x, a)Vu e ■ Vu 



+ (A(x,a £ )-A(x,a))Vu £ -Vu (4.9) 

v p/2 

and by combining the Holder continuity of A(x,t) with Lemma \4- 1\ 

(a E - a)Vu e • Vu < C(J")w(e) || Vw e ||i2 (n) || Vu ||l2(o) 



< CC{F)w{e) \\4>fr , (4.10) 
ffe 2 o(r) 

where C(J-) is a positive constant depending on the constant of regularity T for A(x, t) and C is a positive 
constant which docs not depend on e. For any real numbers p, q with - + | = 1 

a e (x) -A(i,a(i)))Vu 6 • Vu <C(F)\\ a e - a || L p(o) 



• II Vm £ ||i»(n\uj /2 ) II Vu lli«(ny^ /2 ) ■ (4.11) 
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By Meyers' inequality [M] we have that there exists q > 2 such that 

II ||z« ( n\ir„)< C || V« || L 2 (n) (4.12) 
and the same holds for u £ and combining (|4. 10(1 - (14. we obtain 

II A A(.,a(-)) - A A(-,a e (-)) II* = , SU P ( ( A A(.,a(-)) ~ A A(- ,a e (•)) ) ^' ^ 

^effi(r), ||0|| i =1 
«c 2 o(r) 

< C(w(e)+ || a e — a ||zi>(n) ), 

where C is a positive constant independent from e. The above inequality holds for any e, < e < p/2, which 
concludes the proof. ■ 

Proof of Theorem \2.SX Let x S L p and take < s < p/2. We can split the quantity \A(x, a(x))—A(x, b(x))\ 
as follows 

\A(x,a(x)) - A(x,b(x))\ < \A(x, a{x)) - A(x, a E (x))\ 

+ \A{x,a E [x))-A(x,b E (x))\ 
+ \A(x,b £ (x)) - A(x,b(x))\ 

and by the Holder continuity of A(x,t), Lemma T4. II and Theorem 12.21 

\A(x,a(x))-A(x,b(x))\<C(2cj(e)+ || A£ - A£ ||, ), (4.13) 

where C is a positive constant which does not depend on e. By letting e — > + and Lemma T4. 21 we obtain 
the desired estimate. ■ 

Proofs of Theorems 12.41 [2"31 12.61 follow the same line of proofs of Theorems 2.2, 2.3, 2.4 of |AG| by replacing 
the singular solutions of Theorem 13.11 with the singular solutions with compact support in L obtained in 
Theorem 13.51 

4.2 The N-D map. 

The proof of Theorem 12.71 shall be based on the following construction of singular solutions suited for the 
(I2.28|) with local data. The following is well known. 

THEOREM 4.3. Let SI and V be as in Theorem fOl For any r, < r < r , set z T = x° + TV, for 
some x° 6 L p and p, < p < po. If L is the operator of (|3.1[) . with Holder continuous coefficients matrix 
o~ = {o~ij}i t j=i... n , with exponent < f3 < 1, the Neumann's function N a for the boundary value problem 
associated to the operator (|3.ip 

( L N a (x, z T ) = —5{x — Zr), in Q p 
| crVNvix, z T ) ■ v = jg^, on dttp 

has the form 

2-71 

N a (x, z T ) = C n (det(cr(z r ))) _1/2 ( ( T- 1 (z T )(a; - z r ) • (x - z T )j 2 + R{x, z T ), (4.14) 

where C„ is a suitable dimensional constant and the remainder R(x, z T ) satisfies 

\R(x, z T )\ + \x- z T \ \V x R(x, z T )\ < C\x - z T \ 2 - n+a , (4.15) 

for every x £ Q p , \x — z T \ < ro> where C = C(E) is a positive constant depending on E, ro is a positive 
number which depends only on the geometry of Jl and < a < [3. Moreover 

II N a (; z T )\\ m{n) < Cr (2 - n ^ 2 , for any < r < r°, (4.16) 

where C > is a constant which only depends on diam(Q,), X, L, r, h and r°. 
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Proof. See the proof of Corollarv l3.4l and [Mil Chapter 1]. ■ 

THEOREM 4.4. Let Q and T be as in Theorem 1 2. 2\ For any r, < t < tq, set z T — x° + tv, for some 
x° G T p and p, < p < po. If a is the matrix with entries {o~ij}i t j=i... n in (|3.ip and S is an open portion 
of dflp \ dfl with positive distance from dfl, there exists u £ Hl oc (£l p \ z T ) solution to 

Lu = —5(x — z T ), in fl p 
cVu ■ v = 0, on dftp 

aVu-v — — T^r, on S. 

Moreover 

II u(; z T )\\ H1(n) < Ct^' 2 + B, for any < r < r°, (4.17) 
where C, B > are constants which only depend on diam(fl), X, L, r, h and r°. 



Proof. Let N a (-, z T ) be the Neumann function for £1 



div (<t ViVo- (x, z T )j — ~S(x — z T ), in fl p 
aVN a (x, z T ) ■ v = — | , on dft p 

and S be an open portion of dfl p \ dfl with positive distance from dfl. Set 

u(x) = N a (x, z T ) + w{x), for any a; 6 Q p , 
where w € i? 1 (f2 p ) is the solution to 



< 



div (erVw)) =0, in fi^ 

a Ww ■ v = — jQ^-r , on dfl p \ S 



aVwu = -j w 0^, on 5. 



is a solution of the given boundary value problem and by Caccioppoli inequality it also satisfies (|4.17l) . 



Proof of Theorem \2. 7[ It suffices to follow the arguments of the proof of Theorem 12.31 by simply replacing 
the appropriate singular solutions. H 



Remark 4.5. The argument introduced in Theorem \4-4\ o-iso enables to construct singular solutions of the 
type of those introduced in Theorem \3.5\ which however satisfy the zero Neumann condition on dfl \ T. By 
means of such singular solutions it is rather obvious how the proofs of the remaining Theorems \2.2l \2.5\ and 
Corollary 1 2. 61 can be adapted when the local Dirichlet-to-Neumann map is replaced by the local Neumann-to- 
Dirichlet. 
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